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Abstract. We obtain uniform asymptotics for polynomials orthogonal on a fixed and vary- 
ing arc of the unit circle with a positive analytic weight function. We also complete the proof 
of the large s asymptotic expansion for the Fredholm determinant with the kernel sin zj (ttz) 
on the interval [0, s], verifying a conjecture of Dyson for the constant term in the expansion. 
In the Gaussian Unitary Ensemble of random matrices, this determinant describes the prob- 
ability for an interval of length s in the bulk scaling limit to be free from the eigenvalues. 



1 Introduction 

One problem in the random matrix theory is estimation of the probability for a given interval 
to be free from the eigenvalues. In the Gaussian Unitary Ensemble this probability for any 
interval of length 2s in the bulk scaling limit is equal to the following Fredholm determinant: 



The probability of a gap in the spectrum of random matrices from orthogonal and 
symplectic ensembles is also expressed in terms of A(s) (see [Tl 131 ITS]). 



A(s) = det [I — K] 



(1) 



where K is the integral operator on L 2 (0, 2s) given by 
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It was shown by Jimbo, Miwa, Mori, and Sato .2] that (d/ds) In A(s) satisfies a modified 
Painleve V equation (for simpler proofs of this see 011]). 

An interesting question is calculation of the asymptotics of A(s) for large s (the small 
s series are easy to obtain). The first two terms in the expansion of lnA(s) were found 
by des Cloizeaux and Mehta jHj who used a connection with the spheroidal functions. The 
full asymptotic expansion was obtained by Dyson [Bj with the help of the inverse scattering 
techniques for Schrodinger operators. These calculations were partly conjectural. A rigorous 
derivation of the main term was given by Widom [2j using continuous analogues of orthogonal 
polynomials. Finally, Deift, Its, and Zhou 0] (see also that work for a more extensive 
introduction) calculated, as a particular case of a more general result, the full asymptotics 
of the derivative (d/ds) In A(s) using techniques of matrix Riemann-Hilbert problems. This 
settled the question up to the constant term in the expansion of In A(s). The first 3 terms 
in the Dyson expansion are as follows: 



where the constant term cq = (1/12) In 2 + 3£'(— 1), and C( x ) is the derivative of Riemann's 
zeta function. Thus, justification of cq here remained the only problem and it is solved in 
the present paper. Actually, we obtain the first 3 terms: 1 

Theorem 1 The large s asymptotics o/lnA(s) are given by HJ). 

The proof is based on a formula by Deift 9J which connects the determinants of two 
Toeplitz matrices, a formula by Widom jH] for asymptotics of Toeplitz determinants on a 
circular arc, and on asymptotics for orthogonal polynomials on a circular arc which are 
computed here. 

Let f a (9) be a weight function on an arc a < 6 < 2ix — a, z — e te , < a < 7r of the 
unit circle \z\ = 1, and <p n (z,a) = XnZ n + ■ ■ •, n = 0,1,... the corresponding system of 
orthonormal polynomials: 



Such polynomials in the case of the circle (a=0) were first studied by Szego (see [TU] ) 
who, in particular, found several important asymptotics for them as n — > oo. Afterwards, 
asymptotic analysis of such polynomials was carried out by many authors. Specifically 
for the case of an arc (whose study was initiated by Akhieser [H]), see [121 EH EH EE] 
and references therein. However, the full asymptotic expansion at all points z G C for 

1 As this paper was being prepared for publication, an announcement by T. Ehrhardt claiming the same 
result as Theorem 1 (by a different method) was posted on the internet. 

A third solution to the problem by a Riemann-Hilbert approach (related to the present one) is in prepa- 
ration by P. Deift, A. Its, and X. Zhou. 




s 



oo 



(2) 




k, m = 0, 1, 



(3) 
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a wide class of weights became a feasible task only after recent development of Riemann- 
Hilbert problem methods. It was observed by Fokas, Its, and Kitaev ^H] that orthogonal 
polynomials satisfy certain matrix Riemann-Hilbert problems. An efficient method for their 
asymptotic solution (steepest descent techniques) was developed by Deift and Zhou [13 HE] 
and applied for analysis of polynomials orthogonal on the real axis in PHJ EDj (see also 
[2*T] for a different Riemann-Hilbert approach) and on the unit circle in [22] • The case of 
polynomials orthogonal on [—1, 1] (especially relevant for the present work) was considered 
by Kuijlaars, McLauphlin, Van Assche, and Vanlessen |2H1 El] who found full asymptotics 
at all points in the case of a positive analytic weight on [—1, 1] with power-type singularities 
at the end-points. 

In the present paper, we shall give a procedure to obtain full asymptotics for all z for 
polynomials <p n (z, a) and their leading coefficients as n — > oo in the case of a positive analytic 
weight f a {d)- The argument will be similar to that of [23] • We consider 2s /n < a < tt, n > s, 
s — ► oo, which includes both the cases of a fixed arc (a is independent of n) and a varying 
arc. The asymptotics for <j) n (z,a) we obtain are in the inverse powers of nsin(a/2). The 
remainder after k terms is uniform in a. The general solution is given by equations (j61fJM|) . 
The first 2 asymptotic terms for any z can be easily written using ()58I59|) . An example is 
given in ()65I66|) . 

Our solution can be generalized to the following cases: (1) the weight f a (6) has power- 
type singularities at the end-points of the arc (this can be done following [23]); (2) the weight 
f a {9) is not analytic but only smooth enough and positive (one can approximate it then by 
its Fourier series). 

If the weight is symmetric f a (6) = f a (2iT — 8), there exist relatively simple formulas of 
Szego type ^Ql EH1 HEj connecting polynomials on an arc with those on an interval. In this 
case and for a fixed arc one could try to obtain our results from those of [231 121] ■ The present 
argument, however, is more direct. 

For the proof of Theorem 1, we shall only need asymptotics of polynomials with the 
weight f a {9) = 1 arid only at the point e ia . What we need is summarized in the following 
theorem proved in Section 2 (after the argument in the general case is given): 



Theorem 2 Let < a < tc, 7 = cos(a/2) ; 

-ia/2 e ~ia/2 



+ 3.2V ; " 2H K h 

r 2 + = ^—(16-9e ia + 43e~ ta -2e~ 2ia ), r 2 _ = \(-6 + 7e ta - 17e~ ia ), (4) 
3 • 2 y 2 J 

1 / e ia l 2 2/ / \ 1 + 2 cos a . . . 

r_ = cos (a/2), r— , p = n sm(a/2), e > 0. 
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Let f a (6) = 1. Then the polynomial <f> n (z,a) admits an asymptotic expansion for large p in 
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the inverse powers of p. We have for z = e l 



(e ia ,a) = X nl n e ia{n/2 ~ 1/4) ^np~ 



12 3 
y 

1 + — + — + — + \ i 

P P P \P 



(5) 



where \ n is the leading coefficient of <f) n (z, a) = XnZ n + • • ■ for which we have 



Xl-i = T 2n+l 



1 5 ^ ( 1 



An ' 2 5 n 2 



TV 



(6) 



andr^L is a bounded function of a. The derivative of the polynomial <p' n (z, a) = (d/dz)(p n (z,a) 
at e ia can be written as 



n 



ia(n/2-5/4) \' Kl P \- 2 , ia/2 



2 sin a 



[ip 2 + e ia/2 p + r+ 



( r i + r ) + r i e V2 p + ri') + ^(r!zp 2 + r 2 e^ 2 p) + y + O v , 



(7) 



There exists Sq > stzc/i t/iat a// t/ie remainder terms are valid and uniform in a, s, and n 
for a G [2s/n, ir — e], s > so, and n > s. 

Remark The uniformity of the remainders here is crucial for the proof of Theorem 1 . 

After constructing asymptotics for polynomials in Section 2, we give a proof of Theorem 

1 in Section 3. Note that the present method could also be used to obtain the full asymptotic 
expansion of A(s). 

2 Asymptotics of polynomials on an arc 



In the present section we construct asymptotics for polynomials <fi n {z) = XnZ n + ■ • • orthonor- 
mal with a weight f(z) = f a {6) on an arc a < 9 < 2n — a, z = e t9 , for 2s /n < a < n — e, 
e > 0, n > s, s —>■ oo. This includes both the cases of a fixed arc < a < it and the varying 
arc a = 2s/n. The function f(z) is assumed positive and analytic on the arc for a fixed arc 
case, and on the whole circle in the general case. In the general case, we obtain a-uniform 
asymptotics in the inverse powers of nsin(a/2). 

Consider the following 2x2 matrix 



\~ x 6 (z) y- 1 f MUliO^i 

Y( z ) — | Vn^) Xn Je t-z 2ni£ n 



ffU(g) f(Odg I ' 



where £ is the arc a < 9 < 2tt — a of the unit circle traversed in the direction from 2tt — a 
to a, and <p* n (z) = z n (p n (l/z). As is easy to verify (see [22123), Y(z) is the unique solution 
of the following Riemann-Hilbert problem: 
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(a) Y{z) is analytic for z G C \ E. 



(b) For 9 G (a, 2n — a), Y has continuous boundary values Y + (x) as z approaches x = e ld 
from the outside of the circle, and YL(x), from the inside. They are related by the 
jump condition 



YAx) = Y(x) 



1 x n f(x) 
1 



X 



9 G (a, In - a) 



(9) 



Y(z) has the following asymptotic behavior at infinity: 



Y(z) = [l + Ol- 



z" 

: 





—n 



as z — > oo. 



(10) 



(d) Near the end-points of the arc e ±ia 

Y(z) = O 



\z — 




e ±ia \\ 


\z — 




e ±ia \ J 



(11) 



as z 



e ±ia , z G C \ E. 



The function 



= J - + 1 + y/(z-e ia )(z- 



7 = cos(a/2), 



(12) 



which conformally maps the outside of the arc E into the outside of the unit circle, will have 
an important role in what follows (cf. (6.18) of j3]). Here we take the branch of the square 
root which is positive for positive arguments. Note that the boundary values of if)(z), ip+(x) 
as z approaches x G E from the outside of the circle, and from the inside, are related 

as: 

ip + (x)ip-(x) = x. (13) 

Let fi(z) be defined by the equation (J 12)) but with the minus sign in front of the square 
root. Then /i(z) is the mapping of the outside of the arc into the inside of the unit circle. 
Hence \(J>(z)\ < 1, whereas \ip(z)\ > 1 for z G C \ E. Therefore we have that first, for \z\ < 1 
\z/ip(z) 2 \ < 1, and second, for \z\ > 1 



z 




Zfi(z) 2 




z\i[z) 2 




^z) 2 


tfj{z) 2 




ip(z) 2 fi(z) 2 




z 2 




z 



< 1. 



Thus, 



i)(z) 2 

This inequality will be useful later on. 



< 1 



for 



(14) 
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Figure 1: The contour for Riemann-Hilbert problems. 



We now replace the original Riemann-Hilbert problem with an equivalent one which is 
normalized to unity at infinity and has oscillating elements of the jump matrix. Namely, set 

T{z)=r wn Y{z)^{z)- n ^ , <x 3 =(j (15) 

Then, as is easy to verify, T(z) satisfies the same problem as Y(z) (JHJ) but with the changed 
conditions (b) and (c): 



(b) 



(c) 



T + (x) = T_(s) ( xn ^~ 2n xn J}^_ 2n ) , x G E, (16) 



Following the idea of the steepest descent method of Deift and Zhou, we now replace the 
Riemann-Hilbert problem for T(z) with an equivalent one on a system of 3 contours where 
some of the jump matrix elements are exponentially small. Divide the complex plane into 3 
regions as shown on Figure 1 (the contours S 13 lie sufficiently close to S 2 = X for f(z) to 
remain nonzero and analytic in regions 1 and 2) and define the matrix-valued function S(z) 
by the formulas: 

1) in region 1 

S(z) = T(z) ( \ n °), (18) 
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2) in region 2 

S(z) = T(z) ( I ?Y (19) 

3) in region 3 

S(z)=T(z). (20) 

The condition (c) in the problem for S(z) is the same as for T(z), the conditions (a), 
(b) and (d) are different. Namely, 

(a,b) S(z) is analytic in C\(EiUE 2 U£3) with the following jump conditions on the contours: 

S+(x) = S-(x)( j °\ xGEiUEs, (21) 

S4*) = ^)(_ /( ° x) -i /( Q X) ), *eS 2 = E, (22) 



c as z — > oo 



S(z) = I + 0[-\, (23) 



(d) near the end-points of the arc 



log|z — e I log 1 2 — e 
as ^ -> e ±ia , 2GC\(SiUS 2 U E 3 ). 

Recalling Q14j) . we see that, for n large, for x outside some neighborhoods of the endpoints 
of the arc, the jump matrix on Si U E 3 is uniformly exponentially close to the identity. We 
therefore approximate the function S(z) with parametrices inside the mentioned neighbor- 
hoods and in the outside region where we neglect the jumps on EiUE 3 . The Riemann-Hilbert 
problems for the parametrices can be solved and the solution closely resemble that in case 
of polynomials on an interval [23j. 

The parametrix for the outside region is defined as the solution of the following Riemann- 
Hilbert problem: 

(a) N(z) is analytic for zgC\Ee E 2 , 

(b) with the jump condition on E 

N + (x) = N.(x) (flyi /( X) ) ; (25) 
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as z 



oo. 



(c) and the following behavior at infinity 

N(z) = 1 + 

The solution N(z) is found in the same way as in j2B|- Consider the Szego function: 
Viz) 



(26) 



cxp 



y/(z - e ia )(z - e~ ia ) 



2m 



ln/(0 



This function is analytic outside the arc, and its boundary values satisfy 

D + 0)D_0) = f(x), xeE. 

Denote 

Poo = lim T>(z) = exp 

z— >oo 



1 



Then the solution of the above Riemann-Hilbert problem is as follows: 



(27) 

(28) 
(29) 



N(z) = -{V c 



a + a 



-zya — a 
a + a~ l 



-CT3 



1/4 



(30) 



z(a — a x ) 

where the value of the root satisfies the condition a(z) — > 1 as z — > oo. Note that det N(z) = 
1, which allows, in particular, to write a simple expression for the inverse A^(^) -1 . 

Now consider a S- neighborhood Us of the point e ia , small enough so that f(z) is analytic 
and nonzero there. The jump matrices on Ei, S3 are not close to the identity in this region, 
so we need to construct a separate local parametrix. We look for a matrix-valued function 
P(z) which is analytic in Us, satisfies the same jump relations on (Ei U E 2 U E 3 ) fl Us as 
S(z), has the same behavior at z = e ia as S(z), and matches N(z) at the boundary: 

P(z)N(z)- 1 = 1 + ( , ) ) , zedU s \ (S x U E 2 U S 3 ), 

\nsm{a/ 1) J 

where p = nsin(a/2) — > 00. We seek P(z) in the form 

'v(z)^ 



(31) 



P z) = E z)P z) 



-<r 3 /2 



(32) 



where E(z) is invertible and analytic in a neighborhood of Us- The function E(z) does 
not affect jump relations and will be chosen later so that P(z) satisfies (|31|). Using the 
boundary- value property (fT3*j). we obtain (cf. 



P(x)+ = P(x). 
p(x)+ = P{x). 



1 
1 1 

1 
-1 



x E (Ei u s 3 ) n Us, 
x e s 2 n u s . 



(33) 
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Consider the function uj{z) defined by the equation: 

Using Q13jl . we have for the boundary values of to on the arc: 



(34) 



^M + = ln^vZ = - v / ^M_, (35) 



therefore u)(z) is analytic in Us- For z near e ta , we obtain uniformly for all a 



*{z) _ 1+ / ism( a /2) e _ w 2 / 1 + _^^ +(eW2cQs(a/2rl _ 1)e - t ^ | (1 + 0(tt))j 



y'i I y cos(a/2) y 2zsina 2 



(36) 

If also |w| < 2 sin a, we get for the function u)(z) the following (nonuniform) expansion at 
z = e ia - 

,s «sin(a/2) ■ ( 1 - 2e~ ia - 2e~ 2ia ,,\ , , 

u;(z) = r4rv e M 1 — u + 0{u 2 ) . (37) 



cos(a/2) \ 6isina 

Denote 

P(z) = Q((), C = n 2 u(z)- (38) 

Now we reached a crucial moment. The circle dUs is transformed in the ( variable into 
a curve dUg whose minimal distance from zero is n 2 mino<t<27r \uj(e ta + 5e lt )\. In order to 
construct a solution, we need this distance to be large. This is so for large n if S and a are 
independent of n (see (137)1 ). In the general case of 2s /n < a < tc — e, there exists some small 
a (depending on f(z), e) such that we can assume 



^ _ f sin(o!o/2), for a < a < n — s 
1 sin(o;/2), for 2s/n < a < a 



(39) 



Putting u = 5e lt in ()36|) and choosing ao sufficiently small, we obtain after simple analysis: 

n 2 min \uo(e ia + 5e u )\ > C(nsm(a/2)) 2 (40) 

for some constant C which is larger than zero. We see that for large s this distance remains 
uniformly large (not less than of order s 2 ) for any a G [2s /n, n — e], where n > s. 

In Us the image Si,2,3 of the cuts can be considered as 3 direct lines emanating from 
zero. (The image of S is a line, and the exact form of Si and S 3 can be chosen at will.) 
The analytic matrix-valued function Q(() which satisfies the jump conditions (|33|1 on £1,2,3 
and singularity conditions (|24|) at ( = was constructed in [23] ■ Namely, we have in Us 
(the regions 1, 2, and 3 correspond to the (^-variable images of the regions in Figure 1) the 
following expressions in terms of modified Bessel and Hankel functions (see, e.g., [2*7]): 
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1) region 1 



(1)^-^/2^1/2 



Q(0 



2) region 2 



2 UC^KM (e-/ 2 C 



-i7r/2/-l/2\ 



7T 



C 1/2 f^ 2) )'(e-/V /2 ) 



(41) 



Q(C) 



1 

2 



C 1/2 (#£°Y (e^ 2 C 1/2 ) vrC 1/2 (tf^Y (e-/ 2 C 1/2 ) 



3) region 3 



0(0 



/o(C 1/2 ) ^o(C 1/2 ) 

«C 1/2 ^(C 1/2 ) -C 1/2 ^(C 1/2 ) 



(42) 



(43) 



where — 7r < arg(^) < 7r. 

We now have to choose E(z) so that the matching condition (|31j) is satisfied. For that 
we can use the first term in the asymptotic expansion of Bessel and Hankel functions for 
large (. The expansion for Q(() is the same in all the 3 regions. We can write down an 
arbitrary number of terms. Below we shall make use only of the first three. We have: 



1 + 



-2i 



-4z 



+ o(c 3/2 ) 



(44) 



2 7 C \ 4i 1 

uniformly on the boundary of Ug. We now define E(z) as follows: 



E(z) = -LN(z)f( z y^( \ I 



(45) 



We verify exactly as in j2H| that it is an analytic function in Ug. Using (J32~|) . (pMJ) . and ((301) 

(and estimating N(z) as below) we see that the matching condition (j31)l is now satisfied. 
Thus we have 



P(z) = E(z)Q(n 2 u(z)) 



/(*) 



-03/2 



(46) 



Solution in the neighborhood Ug of e^ 2n a ^ is similar (but note the reversed direction of the 
contours). We have there: 



P(z) = E(z)a 3 Q(n 2 u(z))a 3 
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/(*) 



1 7 =N(z)f(zr/ 2 () j)(7rnV^)) CT3/2 - 



(47) 



We are now ready for the last transformation of the Riemann-Hilbert problem. Let 



R(z) = S(z)N-\z), z e C \ {U s U Us U E ll2)3 ), 

R(z) = S(z)P- 1 (z), zeU 5 \ S lj2i3 , (48) 

R(z) = S(z)p- 1 (z), zeU s \E 1A3 . 

It is easy to see that this function has jumps only on dU$, dUs, and parts of Ei, and E 3 lying 
outside of the neighborhoods Ug, Us (we denote these parts S^). Namely, 



R + {x) = H.(x)N-(x) ( }n ? j N-{x)~\ x E E° ? 

R+(x) = R^(x)P(x)N(x)-\ xedU s , 
R+(x) = R^(x)P(x)N(x)-\ xedU & . 



out 
3 > 



(49) 



As is easy to verify, the matrix elements of N(x) and N(z) 1 remain bounded for \x — e ±ia \ > 
5 > sin (s/n) and all a. For example, 



Ia 4 | 



l+87sint+167 2 )" 1/2 < (I-47)- 1 for z = e M +sin(s/n)e !i 6 9[/ sin(s/n) . (50) 



(For evaluation of T>(z) near e ±ia see a similar calculation in Lemma 6.4 of |23j.) Therefore 
the jump matrix on E*^ can be uniformly (both in z and a) estimated by ()34|4()j) as I + 
0(exp (— Cinsin(a/2))), where C\ is a positive constant. Note that in the case of a fixed 
arc, it is sufficient to use (fH|). 

The jump matrices on dUs U dUs admit an asymptotic expansion in powers of 1/vT 
(which turns into expansion in powers of 1/n for a fixed arc or 1/s for a varying arc case). 
First, flSD, (B3), and flU yield 



P^jVOzp 1 = / + N{z)f{z) a3 ' 2 



-1 -2i \ 3/1 -4z . 

.8v? V" 2< lj 2'ci« 1 ) + (51) 

o(r 3/2 )] /(zy^Niz)- 1 = i + a 1 + a 2 + o(r 3/2 ), * e at/*, 



where Ai and A 2 denote the terms with y/C, and (, respectively (the remainder term will be 
justified below). We write them down explicitly for the case of f(z) = 1 needed later on: 

A 1 f"(3« 2 -a- 2 ) -*(3a 2 + a" 2 ) \ A __ 3 fl ( . 9] 

1 ~ 2 4 ^ V -^( 3fl2 + a_2 ) 3fl2 " / ' 2 7 C V 4i 1 / ' 

where / = 1. The functions Ai(z) and A 2 (z) for 2 e dC/a are given by the same expressions 
with a exchanged with 1/a and % replaced by —i. 

Using the expansion for Bessel functions, we can write a general term A 3 - in (|51|) which 
is of order \jQl 2 . Indeed, apart from the prefactor with C^ 2 , the matrix elements of Aj are 
obviously 0(1) as ( — > 00 for a fixed arc. It is also true in the general case, because of the 
remark after equation (|49|). Now it is clear that (|51|> is an asymptotic expansion in (. 
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Since for z G dUg by (|5?l47Ij) 

= n^AT) = 0(n sin(a/2)) = { » ^ (53) 

the component Aj on this boundary is of order l/n J (fixed arc) and 1/V (varying arc), and 
the remainder term in ()51|) is uniform for all a between a fixed positive value and 2s /n, 
all n > s, and all s larger than some sq. (As we shall see this uniformity persists for the 
remainder in the asymptotics of our polynomials.) We show as in j2S] that Aj(z) is an 
analytic function in Ug with a pole at e ia of order less than or equal to [(J + l)/2]. The same 
reasoning also holds for the neighborhood Ug. We shall denote the components of the jump 
matrix there by the same symbols Aj(z). 

If R(z) is known, we can trace the sequence Y i— ► T i— > S i— > i? backwards, and obtain 
an expression for the polynomials. 

We look for -R(z) asymptotically in the form R(z) ~ Rq(z) + Ri(z) + Rj2{z) + •••■, where 
Rj{z) is of the same order as Aj (in our case, of order (nsin(a/2))~ J ). For more discussion 
and justification of this expansion see ^H1I23I- More precisely, it can be shown as in Theorems 
7.8-7.10 of PH that for any k > 1 

fc-i 

R(z) = I + J2 R i( z ) + 0((nsm(a/2))- k ) (54) 

3=1 

uniformly for all z if the arc is fixed, and for z outside a neighborhood of z = 1 if the arc 
is varying (a = 2s/n). The proof of this expansion in the general case 2s/n < a < tc for z 
close to 1 requires a special argument since the contours then are close to z and can not be 
trivially deformed. Such a proof will be given later on for the case of / = 1, z G Ug, we need 
below. The argument in the general case is similar. Moreover, it follows directly from the 
proofs that the remainder term in (|54j) has the same uniformity property as that in (J51j) . 

Substituting this asymptotic expansion into and collecting the terms of the same 
order, we obtain: 

i?o+(s) + Ri+(x) + (Ro-(x) + Rx-(x) + ■ ■ •)(/ + Ai(a;) + ■•■), xedUgU dUg. 

R +(x) = R -(x) Rq(z) = I, 

Ri + (x) — Ri-(x) = Ai(x), 

R 2+ (x) - R 2 -(x) = R 1 ^(x)A 1 (x) + A 2 (x), 

k 

R k+ {x) - R k -(x) = ) j R k -. j -(x)A j (x), k = 1,2,... 

3=1 

(55) 

The main term in the asymptotics of polynomials is given therefore by the parametrices 
at the appropriate points z. The expressions for Rk{z) follow from the Sokhotsky-Plemelj 
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formulas: 

2?™ JdUgUdUs x — z 2ni Jd UsU du s x — z 

(56) 

Note that the contours are traversed in the negative direction (see Figure 1). 

Following we can also obtain the expressions for Rj{z) in a different way. As 
mentioned, Ai(z) is analytic in Us U Us with simple poles at the end-points of the arc. Thus, 

Mz) = ^-- + 0(l), asz^e™, A 1 (z) = + 0(1), as e"*, (57) 

where the constant matrices and B' 1 ' are obtained by expanding various functions in 
()51|) at z = e lOL and z = e~ ia . It is easy to verify directly that the Riemann-Hilbert problem 
for R\(z) has the solution: 



A W . B (i) 



^ + for *e<C\ (17,1^) 



^ + y _ g - <a - Ai(z), for z e Us U U s . 



z—e l " z—e 



Expanding the functions in (|5ip. we obtain: 

AW = C ° S(a/2) f .I 2 "^W, 5 (1) =^), (59) 

where M means complex conjugation applied to every matrix element of M. The general 
term Rk(z) is obtained similarly provided we have the expressions for Rj(z), j — 1, 2, . . . , k — 
1. Since Rj(z) are analytic in UsUUs and Aj(z) have poles at e ±ia of order at most [(j + l)/2], 
we see that 

k aW A (k) Ah) 

R k -i Jz)A j (z) = 7 ^^ + 7 Ez V^ + "-- + — — + 0(1), asz^e", (60) 

where p = [(k + l)/2], and similar expressions hold with matrices A replaced with some 
matrices B in a neighborhood of e~ ta . Then 



EJ=i ( + (7-^)7 ) , for * G C \ (U s U f/ 5 ) 

E? =1 ( + (^f) - E -=i i2fc- i ,-(z)A i («) J for z & Us U U$ 



(61) 

Recalling the definitions of i£, S, and T, we finally get 

Y(z) =^R(z)M(z)iP(z) nr7 ^ R(z) ~ I + R^z) + R 2 (z) + ■ ■ ■ , (62) 
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where for z restricted to region 3 M(z) = N(z), P(z), or P(z) if z G C \ (Us U Us), Us, or 
f/,5, respectively (the expressions in regions 1 and 2 can also be readily written). Therefore, 
by©, 



Y 11 (z) 



y n ^ n (z)[R n (z)M n (z) + R l2 (z)M 21 (z)\. (63) 
Furthermore, 

xLi = >2i (0) = 7~ 2n [#2i(0)iVn(0) + i? 22 (0)iV 21 (0)]. (64) 

Since we know the expressions for M(z) and can obtain R(z) with an arbitrary precision, 
the last 2 equations give an implicit solution for asymptotics of polynomials <p n (z) and their 
leading coefficients. These are the asymptotic series in the inverse powers of nsin(a/2). The 
error after k terms is 0(n sin(a/2)) _fc_1 and remains uniform for all a, s, and n, provided 
a G [2s /n, it — e], n > s, s > so. 

As an example, we give below the first two terms in the asymptotics of <p n (z) valid for 
z outside a fixed arbitrary small e-neighborhood of the arc E: 



7 > n (z) Vs 



V(z) I 



a(z) + a(z) 



8n 



a(z)e 



ia/2 



z — e 1 ' 



+ 



a(z) 



-lg-ta/2 



X 



1 



-2n+l 



n-1 



1 

An 



(65) 
(66) 



V(0)V x 

where all the quantities are defined as above (see (fT2l I2T1 12UI IHOjl ). 

We now give an explicit solution for the first 3 terms in the asymptotics of (j) n (z) and 
its first derivative at an end-point of the arc for the weight f(z) = 1. 

Recall that we still have to prove the remainder term in the expansion (|54jl for z G Us, 
2s/n < a < aQ. Since f(z) = 1 we see that T>(z) = 1. Each matrix element of Ak(z) can 
be written as (a^a 2 + (3ka~ 2 + 7fc)C _fe ^ 2 ) where a&, (3k, and 7^ are independent of a. In a 
neighborhood of e ta , we have the series: 



c 



-k/2 



[n usma) 



3=0 



a 



u 



sin a 



a[z) 



u 



sma 



3=0 



u 



sin a 



(67) 



\u\ < sin(a/2), 



u 



z 



< \cj\ < 1, 



J>0, 



where bj(k,a) are bounded functions of a. Since A&(z) is single-valued in its matrix 
elements do not contain terms with y/u. Hence, ack, Pk = if k is even, and 7^ = if k is 
odd. We have in the same neighborhood: 



A,(z) 



1 



nsmm 



-[(fc+l)/2] 



sin J a 



(68) 



Now using (|50|l and the second formula in (jpTj) . it is easy to show by induction that 

1 a. , U 3 



Ruiz) 



(nsin a) 



a: 
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sm J a 



(69) 



The matrix elements of Cj and Cj are bounded functions of a. Considering the remainder 
term in the asymptotic expansion of Bessel functions, we see that R(z) — I — Y^jZi Rj{ z ) is 
given by the same series (jfiH]! with different matrices Cj(k,a) (but also bounded in a). By 
analyticity, these series and their derivative w.r.t. u converge in Us- This is also true for 
a > a . Thus, 



fc-i 



R( z ) = I + ^R j (z) + 0{(nsm{a/2)y 



^-R( z ) = Y^R 3 (z) + 0{n- k sm(a/2)- k - 1 ), z G U s , 

dz dz 

where the remainder terms are uniform for z G U$, 2s /n < a < tt — e, n > s, s > sq. 

Using (j63l46|45|43|) . we obtain for z in the intersection of Us and region 3 an asymptotic 
equivalence: 



(Rj,n(z) + iR j>12 (z))aI' }) , 



^/o + al' + J2i( R iMz) - iR jjl2 (z))a- 1 I + 



3=1 



where the Bessel functions Iq and Jq = I\ are taken at nyu(z). We now estimate the U R' 
terms and some of their derivatives at the point e lOL . From (J5fij) and (J52j) we get: 



2?rz 



3a 2 dtt 1 a 2 
+ Res 7 



27ri J M 2? sin a x e 8\/C' 

where u = 2 — e* a = ee 2< describes a circle of a small radius e in the positive direction. The 
expansion of £ near e m is given by (|37)1 . and that near e~ m is obtained from it by changing 
the sign of a. Expanding also a 2 (z) near e ±lOL and calculating residues, we obtain 

o p ia/2 _i_ p -ia/2 

Rx n (e ia ) - ii?! 12 (e ia ) = — -- . (72) 

' V ; ' V ; 16msin(a/2) V ; 

Similarly, we calculate 

1 4. p-ia _ Opia 

R 1A1 (e ia )+iRiMn = e-* a/ \ + 6 . 6 ■ (73) 

3 • lorn sm(a/2) 

Now differentiating (J55j) w.r.t. z, we get 

Ax (a: 



2tu J a 



dUsUdUs V x ~ z ) 
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From here we obtain as above 



Rl ( P ™\ it? (r^- 6 i<X/2cOS ( a / 2 ) (7A \ 

lonsm (a/zj 

To estimate R 2 (z) using (jSEJ), we need to know L(z) = Ri(z)Ai(z) + A 2 (z) in the 
neighborhoods of e ±ia . Here Ai(z), A 2 (z) are given by (jo^j) and for we use the second 

formula in ()58j) . Then we obtain in the same way as above: 



# 2 ,ii(0-^2,i 2 (e iQ ) 

#2, 11 (0+^2,12(0 



-6 + 7e* Q - 17e~ ia 
2 9 n 2 sin 2 (a/2) ' 
16 - 9e iQ + 43e" ia - 2e 
3 ■ 2 9 n 2 sin 2 (a/2) 



(75) 



Substituting these expressions into (fTTj). expanding 0(2;), Bessel functions, and u;(z) at 
e lct (see (J21J)), and recalling (fTUj). we obtain A simple calculation of (jo^|) yields ©• 
Taking the derivative of ()71)1 at z = e ta and using the expressions for R(e ia ) and R'(e ia ) we 
complete the proof of Theorem 2. □ 



3 Proof of Theorem 1 

Consider the following weight function on the unit circle: 

J w [0, otherwise. 

and the corresponding orthonormal polynomials a) = ^/c^+^-i^ - ^- • -+^o satisfying 
(j2J). Since f a {6) = f a (2ii — 9), the coefficients of the polynomials (f>k{z, a) are real. 

Associated with f a (8) is an (n + 1) x (n + 1) Toeplitz matrix T n (a) whose matrix elements 
are as follows: 

(T n (a)) jh = -±- j** UO)d0 = ( 1 ^Ju-k)) \ = , J' , j, k = 0, 1, . . . , n. 

2k Jo I n (j-k) ' r K 

Putting a = 2s /n and taking the limit n — > oo, we easily obtain 

A(s) = lim detT n _! f — j . (76) 

n-»oc \ n J 

If a is fixed, the large n asymptotics of detT n (a) were obtained by Widom [Bj (see also [2B] 
for an alternative derivation), namely, 

det T n _!(a) = cos n2 (a/2) (n sin |) 2 1/1 V C ' ( - 1) (1 + o(l)), < a < tt. (77) 
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An idea of Dyson jB] was to put a = 2s /n in these asymptotics. Then, in view of (|76|). one 
formally obtains the first 3 terms of (J2J). However, as noted in j^j, since the remainder term 
in (|77j) is unknown, we cannot say if this expansion is uniform in a and therefore cannot 
justify such a limit. 

Recently, Deift ]Hj found a formula which connects the determinants of two different 
Toeplitz matrices. A proof of it we need below was noticed by Simon. We shall use this 
formula together with equation (J77)l for det T n _i(a) (a fixed) to obtain an expression for 
det T n -i(2s/n). The variant of Deift's formula we shall use is the following: 

±\ndetT n _ 1 (a) = -J2±- [** " -^ fc (e* a)\ 2 d9, n = l,2,... (78) 
(X(y. ^ Ziix 3 ot cicx 

Indeed, it is well-known (e.g., jTUJ ) that the Toeplitz determinant has the following repre- 
sentation in terms of the leading coefficients of <f> n (z, a): 

n-1 

detT„„ 1 («) = JJxfc 2 , n = l,2,... (79) 

fc=0 

Therefore 

Al n detT„_ 1 («) = -2X:^. 
On the other hand, the orthogonality of our polynomials implies 

r l ~ j- 1 1*171 — OL J 

n-1 1 r 2n-a n-1 



» i pZTr—a i t 

- E / W e * e > + b' k _M)e~ l{k - 1)e + •••) + cc. ) d9 = -2J2 

k=o Z7r Ja • - 

and (f75|) is obtained. A corollary of it is the following 



LEMMA 3 Let <f) n (z,a) and T n _i(a) be defined as at the beginning of the section, and 
4>' n {z, a) = (d/dz)(f) n (z, a). Then for any n = 1, 2, . . . 

(1 71 ! 

— lndetT n _ 1 («) = -|0 n (e ia ,a)| 2 --{0„(e~^,a)e l Vn(e ia ,«) + c.c.}. (80) 

da 7T 7T 



Proof From the identity 
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we obtain using ([7S|). 

-lndetT^a) = -- V|0 fc (e^a)| 2 
da 7r L — ' 

k=0 



As is known (e.g., |10j). any system of orthonormal polynomials on the circle satisfies an 
analogue of the Christoffel formula which we write for x and y on the unit circle in the form: 

a ( \in~\ ( KM ( ^n(y) - M x )Mv) 4> n (x){(f>n{ x ) - <f>M) - 0n(^)(0n(^) - <ftM) 

fc=0 



1- xy x (x - y) 



where <fi* n (x) = x n <f) n (x). Letting y — > x along the unit circle and noting that (d/ dx)4>* n (x) 
nx~ n+l d>„(x) — x~ n+2 cb'(x), we obtain 



j~] \<pk{x)\ 2 = x4> n {x)(p' n (x) + -4> n (x)4>' n (x) - n\4> n (x)\ 2 , \x\ = 1, (81) 

k=0 

which, after substituting x = e lOL and recalling that the coefficients of our polynomials are 
real, completes the proof. □ 

Note jU |2H] that the logarithmic derivative (J8Uj) satisfies a r-function version of Painleve 
VI equation. 

Now we want to integrate (jHUjl over a between a± = 2s /n and any fixed «2 < vr for large 
n (for a?2 asymptotics ()77|) are valid). We therefore need to know large n asymptotics of 
4> n {z, ®) and its derivative at the endpoint z = e ia of the orthogonality arc. These are given 
by Theorem 2. Substituting them into (jHUjl . we see that the term ^|0 n (e ia , a)\ 2 cancels at 
once, and purely imaginary terms in the rest of the expression also disappear (in particular, 
the terms with drop out). As a result we have: 

d , , m / s 9 sin(a/2) cos(a;/2) ^ ( 1 \ . „. 

— lndetT n _ x (a = -n 2 \' - . \ 1 ' + O . 2 , 82 

da 2cos(a/2) 8sm(a/2) \nsin (a/2) J 

where the remainder term is uniform for 2s /n < a < tt — e, s > s , n > s, e > 0. Integrating 
this expression over a from 2s /n to a2 and using Widom's asymptotics (fTTjl for det T n _i(a 2 ), 
we get 

si s ( 1 \ 

In det T n _i(2s / n) = n 2 In cos Inn sin — h c + O I — -— I + o(l) (83) 

n 4 n sin(s/n) J 

as n — > oo with the first remainder term turning into 0(1/ s) valid for all s > So- By f!76|) . 
this equation yields (0). □ 

Remark The Riemann-Hilbert problem methods allow us to calculate asymptotics of or- 
thogonal polynomials to arbitrary precision. Because of the integral identity (Theorem 2b 
and equation (4) of (201, see a ^ so Ell) that expresses Toeplitz determinants in terms of or- 
thogonal polynomials, an arbitrary number of terms in asymptotics of Toeplitz determinants 
could also be computed. Thus one should be able to estimate the remainder term in (|77ji 
and therefore give another proof of Theorem 1. 
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